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ABSTRACT 

Using the Galactic A code of Benson et al., we obtain quantitative measurements of 
spheroid-to-disk ratios for a sample of 8839 galaxies observed in the Sloan Digital Sky 
Survey. We carry out extensive tests of this code and of Gim2D, finding that they 
perform similarly in all respects. From the spheroid and disk luminosities, we construct 
luminosity and stellar mass functions for each component and estimate the relative 
luminosity and stellar mass densities of disks and spheroids in the local Universe. 
Assuming a simple one-to-one mapping between between spheroid mass and the mass 
of a central supermassive black hole, we provide the most accurate determination so far 
of the black hole mass function in the local universe. From this, we infer a cosmological 
mass density of black holes of p. = (3.77 ± 0.97) x lO^/iMoMpc"^. We compare our 
results to predictions from current hierarchical models of galaxy formation and these 
are found to fare well in predicting the qualitative trends observed. We find that stars 
in disks contribute 35-51% of the local stellar mass density. 

Key words: galaxies: structure; galaxies: abundances; galaxies: bulges; galaxies: 
luminosity function, mass function; galaxies: statistics 



1 INTRODUCTION 

The distinction between disks and spheroids is one of the 
defining properties of galaxies. Determining the relative im- 
portance of these two basic types of galactic component is 
fundamental to a broad characterization of the galaxy pop- 
ulation. Yet, this is a complicated task which requires not 
only high quality imaging for large samples, but also soft- 
ware capable of decomposing the light from each object into 
a disk and a spheroiqj. The determination of spheroid lu- 
minosities has recently received even more prominence since 
the discovery that perhaps all galaxies harbour a supermas- 
sive black hole at their centre whose mass is proportional to 
the luminosity of the spheroid or bulge (Kormendy & Rich- 
stone 1995; Magorrian et al. 1998; Merritt & Ferrarese 
2001; Marconi k Hunt 2003; Haring & Rix 2004). 

From a theoretical point of view, explaining why most 
of the stars in the Universe end up either in disks or in 
spheroids and understanding the physical processes that re- 
sult in the formation of one or the other of these morpholog- 



^ We will use the term "spheroid" throughout to refer to both 
elliptical galaxies and the bulges of spiral galaxies. 



ical structures is a major challenge. The current theoretical 
framework used to investigate galaxy formation is the cold 
dark matter model (Peebles 1982, Blumenthal et al. 1984, 
Davis et al. 1985), in which galaxies build up hierarchically. 
Within this model, the basic processes thought to be re- 
sponsible for the distinction between disks and spheroids 
were identified over twenty years ago (Fall 1979, Frenk et 
al. 1985): disks result from the collapse of rotating gas cool- 
ing within dark matter halos whereas spheroids result from 
major mergers or disk instabilities (Fall & Efstathiou 1980; 
Barnes k Hernquist 1992; Mo, Mao & White 1998). The 
traditional categories of galaxy morphology, spirals, irregu- 
lars, etc, are too detailed for current theoretical models to 
explain, but the relative luminosities and stellar masses of 
spheroids and disks can readily be predicted (Kauffmann, 
White & Guiderdoni 1993; Baugh et al. 1996a; Baugh et 
al. 1996b; Kauffmann, Chariot & White 1996; Kauffmann 
& Chariot 1998; Somerville & Primack 1999; Hatton et 
al. 2003). Thus, accurate measurements of these quantities, 
for example, as a function of absolute magnitude and in dif- 
ferent environments provides a powerful test of models of 
galaxy formation and evolution. 

An early attempt to determine the relative contribu- 
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tions of spheroids and disks to the luminosity density of 
the Universe was made by Schechter & Dressier (1987). 
They studied a magnitude limited sample of ^ 200 galax- 
ies brighter than V — 16.5, drawn from the catalogue of 
Dressier (1980a), and determined spheroid-to-disk ratios by 
visual inspection. From this, they derived the distribution of 
spheroid-to-disk ratios, as a function of absolute magnitude, 
and found the overall spheroid-to-disk ratio to be higher in 
high density environments (galaxy clusters) than in low den- 
sity environments (the "field"). Schechter & Dressier (1987) 
found that disks appear to contribute roughly twice as much 
as spheroids to the mean luminosity density of the Universe. 
Since a large fraction of the disk light comes from a relatively 
small number of young stars, Schechter & Dressier (1987) 
concluded that the relative contribution of spheroid and disk 
components to the mean stellar mass density of the Universe 
is very nearly equal. 

More recently, Benson et al. (2002) developed a quanti- 
tative method to determine galaxy morphology, specifically 
to estimate spheroid-to-total (S/T) light ratios. This method 
is implemented in the code Galactica (GALaxy Auto- 
mated ComponenT Image Construction Algorithm) . Benson 
et al. (2002) analyzed a magnitude-limited sample of ^ 100 
field galaxies brighter than / = 16.0 and found the luminos- 
ity functions of spheroids and disks to be remarkably simi- 
lar. They provisionally concluded that spheroids and disks 
contribute almost equally to the total stellar mass density 
in the Universe but stressed the significant uncertainties in 
their result arising from the small sample size. A larger sam- 
ple of 1800 galaxies drawn from the Sloan Digital Sky Sur- 
vey (SDSS; Stoughton et al. 2002) was analyzed by Tasca 
& White (2005) using the pubhcly available code Gim2D. 
They found that 54 ±2% of the local cosmic luminosity den- 
sity in both r and i bands comes from disks, 32 di 2% from 
"pure bulge" systems and the remaining 14 ±2% from bulges 
in galaxies with detectable disks. 

In this paper, we perform spheroid/disk decompositions 
from r-band images of a much larger sample of galaxies in 
the SDSS. The structure of the paper is as follows. In 321 
we describe our dataset and how it is processed. In 331 we 
present results from the spheroid/disk decomposition. In 311 
we derive the luminosity functions of disks and spheroids, 
and present stellar mass functions and also the supermassive 
black hole mass function. Finally, in 33 we give our conclu- 
sions. Appendices describe extensive tests of the reliability 
of both Galactica and Gim2D (Appendix IA|) . technicali- 
ties of the fitting process (Appendix [BJ and comparisons of 
our S/T ratio morphologies with more traditional morpho- 
logical measures (Appendix [Cj. A cosmological model with 
Qo = 0.3, Ao = 0.7 is adopted throughout and the Hubble 
constant is defined to be Hq — 100/ikms~^ Mpc~^. 



2 DATA: SLOAN DIGITAL SKY SURVEY 
2.1 Basic properties 

The Sloan Digital Sky Survey (SDSS) is the largest imag- 
ing and spectroscopic survey to date. The SDSS Early Data 
Release (EDR), made publicly available in 2001, consists of 
a 462 square degree area imaged in five passbands {u,g,r,i 
and z) and also covered spectroscopically. The SDSS EDR 



galaxy catalogue is spectroscopically complete down to r = 
17.7 and contains measurements of various galaxy parame- 
ters (Stoughton et al. 2002). The imaging data were taken 
with a dedicated 2.5m telescope in the drift-scan (time- 
delay) integration mode with an effective exposure time of 
54s. The data used in this study are the r-band imaging 
frames with corrections for bias, flat fleld, cosmic ray and 
pixel defects (Lupton et al. 2001). Each imaging frame is a 
2048 X 1489 pixel array with a pixel size of 0.394''. 



2.2 SDSS apparent magnitude limit 

Benson et al. (2002) measured spheroid-to-total (S/T) light 
ratios for the field galaxy sample of Gardner et al. (1996) 
using /-band imaging. The data were originally obtained to 
determine the K-hand luminosity function and Benson et al. 
(2002) showed that they could be used reliably to estimate 
S/T ratios for galaxies brighter than /car = 16.0 with an 
rms accuracy of arms 0.1. Unfortunately, the area covered 
by this sample is rather small (4.4 deg^). 

The SDSS imaging data were obtained with a larger 
telescope but using shorter exposure times than those of 
Gardner et al. (1996). From Monte-Carlo simulations, Ben- 
son et al. (2002) established the signal-to-noise required 
to obtain reliable measurements of S/T using the Galac- 
tica decomposition code. Assuming that, when applied to 
the SDSS data, the code will be reliable to the same overall 
signal- to- noise level, we find that the limiting magnitude re- 
quired for our SDSS sample is /sdss — ^Gar = 0.4. Using the 
mean galaxy colours of Fukugita et al. (1995), the transfor- 
mation between the /sdss and r bands is r — /sdss = 0.9, 
making the total difference equal to r — /oar = 1.3 magni- 
tudes. We therefore select EDR galaxies with r ^ 17.3 — this 
is 1.4 magnitudes fainter than the sample used in a similar 
study by Tasca & White (2006). 



2.3 SDSS data selection and galaxy catalogue 

Galaxies with r ^ 17.3 in the SDSS EDR equatorial strip 
are plotted in Fig. [1] colour-coded according to the SDSS 
run number (94, 125, 752 and 756). The black points repre- 
sent imaging taken in 'poor' seeing conditions (PSFfwhm > 
1.55^^ where PSFfwhm is the full- width at half-maximum 
of the point spread function). This cut on of the seeing is 
used to impose a second galaxy selection criterion since reli- 
able spheroid-to-disk decompositions require that the seeing 
be less than a typical galaxy half-light radius (Beijersbergen 
et al. 1999). 

The final galaxy selection criterion is redshift. To avoid 
contamination of the measured redshift by the local galaxy 
infall velocity, a low redshift cut, z = 0.02, is imposed. Since 
the total SDSS sample begins to tail off at large distance, a 
high redshift cut, z = 0.3, was also imposed. 

The selection leads to a total of 8839 SDSS EDR galax- 
ies. 



2.4 Sample solid angle 

To calculate the solid angle covered by our sample, the 
galaxy coordinates were accumulated in 0.2° bins. All the 



Luminosity and stellar mass functions of disks and spheroids in the SDSS and the SMBH mass function 



V 






Figure 2. The radial Moffat profile {(3 = 4.5) fits (solid line) to 
the SDSS stellar light profiles (points) for stars found at various 
positions within several SDSS frames. The radius is in pixels and 
the Pixel Value are counts. Results are shown for two different 
frames. 



Figure 1. SDSS galaxies that meet our selection criteria. Black 
points correspond to galaxies imaged when the seeing was greater 
than 1.55'^ 



areas which contain at least one galaxy residing in a par- 
ticular bin were summed to give the total solid angle. The 
bin size was chosen such that the derived solid angle was 
insensitive to small changes in the bin size. As an additional 
check, this same bin size was used to reproduce the SA of the 
entire SDSS EDR. For our chosen bin size, the solid angle 
subtended by our sample is 165.5 square degrees. 

2.5 Object detection and astrometry 

Object detection was performed using S Extractor v2.2.2 
(Bertin & Arnouts 1996). The SExtractor world coor- 
dinates of the object centroid positions {x^y) were used to 
identify the catalogued galaxies within the SDSS frames. 
The Galactica code (Appendix IA2.3|) was run on the ex- 
tracted postage stamps whose size was set equal to {2xRp)x 
(2 X Rp), where Rp is the Petrosian radius (Lupton et al. 
2001). This is large enough to contain many background pix- 
els but sufficiently small to ensure a reasonable convergence 
time for the fitting procedure. Prior to decomposition, the 
SExtractor estimate of the local sky background was sub- 
tracted from every postage-stamp to ensure that the back- 
ground level was close to zero (see Appendix lA)) . 

2.6 SDSS point spread function (PSF) 

Before beginning the decomposition procedure, it is neces- 
sary to ensure that the PSF analytic model assumed by the 
Galactica code (Appendix I A2.3|) is a realistic representa- 
tion of the SDSS PSF. To demonstrate that the SDSS stars 
are well represented by the analytic Moffat profile assumed 
by the Galactica code, the iraf task IMEXAMINE was used 
to fit radial Moffat profiles to a sample of stars imaged on 
different SDSS frames and at different positions within ev- 
ery frame. Fig. [2] shows radial fits to stellar light profiles 
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Figure 3. The variation of the PSF across an SDSS frame. The 
dashed line represents the mean value of the seeing for the frame. 
The seeing appears not to vary by more than ±5% from the mean 
value. A similar inspection of other frames shows this to be true 
in general. 



obtained using /3 = 4.5 and demonstrates that a Moffat star 
with this value of /3 is a good analytic representation of the 
SDSS PSF. 

2.6.1 PSF variation 

The Galactica code assumes a starting value for the PSF 
equal to the measured value of the seeing in the SDSS and 
allows the value to fluctuate by ±5% (Appendix IA2.3|) . The 
±5% variation is set from the observed variation of the see- 
ing across a typical SDSS frame, as shown in Fig. [S] which 
demonstrates that for stars imaged at various positions in 
a given SDSS frame, the FWHM does not change by more 
than ±5%. Galactica assumes a circularly symmetric PSF, 
although this may not be precisely true for drift scan obser- 
vations such as those of the SDSS (Berstein & Jarvis 2002). 

The small allowed change in the seeing ensures that 
the Galactica code can find the representative value of 
the seeing at each galaxy position. However, it is important 
to test how consistently the Galactica code recovers the 
'correct' representative PSF for a given galaxy and quantify 
the effect this has on the recovered S/T ratios. The observed 
galaxy properties are expected to vary little between the r 
and i bands but the PSF signatures for these observations 
will be somewhat different. Fig. |4] shows a good correlation 
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Figure 4. (a) Correlation between S/T ratios obtained for the same galaxies in the r and i bands. The good correlation demonstrates 
that the S/T ratios are accurately determined for different representative PSFs and across the apparent magnitude range. The recovered 
S/T ratios show no obvious dependence on galaxy apparent magnitude, indicating that the decompositions are not affected by the 
variation in the signal-to-noise ratio, (b) The difference in the output Galactica seeing for the same set of galaxies observed in the r 
and i bands. The lack of a trend demonstrates that the Galactica code recovers the representative PSF for each galaxy well, without 
biasing the recovered S/T ratios. 



between the S/T ratios obtained for the same set of galaxies 
imaged in the two bands. The code finds consistent S/T 
ratios across a range of apparent magnitudes independently 
of the seeing. Error bars are obtained from 30 Monte Carlo 
realizations of each of the model fits, assuming the noise 
appropriate to each image. 



3 GALACTICA DECOMPOSITIONS AND 
GALAXY MORPHOLOGIES 

3.1 Galactica decomposition outputs 

Fig.[5]demonstrates a typical fit to a galaxy light profile. The 
figure shows the postage-stamp of a real galaxy, a noise- 
free model generated from the best-fit parameters, along 
with the individual model disk and spheroid components. 
In this study, a galaxy is deemed to be sufficiently well rep- 
resented by the model if the per degree of freedom satis- 
fies < 2.0, and if there are no obvious structures left in 
the residual image. Obviously for well- fit data and correctly 
estimate noise should be very close to unity. Allowing 
for larger values of x^ allows galaxies with small (yet sig- 
nificant) departures from our photometric model to be in- 
cluded in our final sample. An example of a well fit galaxy 
is shown in Fig. (6] The cross-hatched areas represent poten- 
tial contamination from overlapping objects as determined 
by the Galactica masking procedure (see Appendix IA2.3|) 
and are excluded from the fitting. The inset in Fig.[6]shows a 
histogram of dP/d(5'/T) — the distribution of the spheroid- 
to-total ratio from 30 Monte Carlo realizations, with the 
vertical dashed line indicating the best-fit S/T value for this 
galaxy. 



S/T= 0.30 ±0.07 

Figure 5. Top: real (left) and model (right) images. Bottom: disk 
(left) and spheroid (right) component fits. The cross-hatched re- 
gions represent potential contamination from overlapping objects 
(or regions where data were unavailable after the image was re- 
centred by Galactica) and are excluded from the fitting. The 
contours indicate the pixel values in ADU/s. 



3.2 Correlations of S/T with other fit parameters 

Understanding the properties of this large statistical sample 
of galaxies is important since it may reveal features which 



Luminosity and stellar mass functions of disks and spheroids in the SDSS and the SMBH mass function 5 



S/T= 0.30 ±0.07 

Figure 6. Real (left) and residual (right) images. The inset shows 
the distribution of S/T ratios from 30 Monte Carlo realizations 
with the vertical dashed line indicating the best-fit S/T. The value 
of is acceptably small and the residual image also shows a good 
fit to the data. 



otherwise would not be discovered in smaller samples such 
as those discussed in Appendix [C] Equally, any unexpected 
correlations between parameters could help discover and re- 
duce possible biases introduced by the fitting routine. 

Histograms of various properties of our SDSS galaxies 
inferred from the Galactica decompositions are shown in 
Fig. [71 These plots reveal the following: 

(i) a large number of highly elliptical spheroids; 

(ii) an excess in the number of galaxies with spheroid po- 
sition angle, ^s, equal to 0° and 180°; and, 

(iii) a non-uniform distribution of the cosine of the disks' 
inclination, cos(i). 

In the remainder of this section, we explore the possible 
origins of these unexpected distributions and their influence 
on the recovered values of the S/T ratio. 

3.2.1 S/T vs. Ellipticity 

Around 15% of galaxies appear to have a highly elliptical 
spheroid component whose ellipticity has reached the im- 
posed upper limit|3 of e — 0.83. A large number of frames 
have been inspected by eye and show that these galaxies 
generally exhibit bar-like structures in the direction of the 
detected highly elongated spheroidal component. In these 
cases, the existence of this extra component, which is not 
part of the fitted model, drives the code to fit small and 
highly elliptical spheroids (Fig. [8|). While in principle it is 
possible to include additional components, such as bars, in 
the photometric model (see, for example, Gadotti & Kauff- 
mann 2007) we have not done so here due both to the fact 
that fitting them would result in prohibitively long times 
to fit each image and that, for poorly resolved galaxies, 
additional components can cause further systematic errors 
such as the one described above. These galaxies are gener- 
ally disk-dominated (with mean S/T ratio of 0.14) and thus 
we expect this shortcoming of the model to introduce only 
a small bias on the overall S/T ratio. (It should be noted, 
however, that this problem may be occurring even in cases 

^ The upper limit for the ellipticity corresponds approximately 
to that of the most elliptical observed galaxies (Lambas, Maddox 
& Loveday 1992). 



S/T= 0.37 

Figure 8. An example of a galaxy with a highly elliptical 
spheroid. This galaxy demonstrates how the central bar-like struc- 
ture in the galaxy results in the detection of a highly elliptical 
spheroid along the same direction. The top row shows the real 
(left) and model (right) images. The middle row shows the disk 
(left) and spheroidal (right) components. The bottom row shows 
the real (left) and residual (right) images. 



where the fitted ellipticity is less than 0.83 if seeing has 
made the bar component appear rounder.) Any bias that 
is introduced would increase this ratio, resulting in a slight 
overestimation of the spheroid luminosity density in 21 

Similar problems of this type (i.e. fitting of additional 
photometric components of galaxies such as bars or isopho- 
tal twists by a component of the photometric model) have 
been noted and discussed by Tasca & White (2006) and 
Simard et al. (2002). In such cases, the S/T ratio will be 
incorrectly estimated. We return to this problem in 21 



3.2.2 S/T vs. spheroid position angle 

Many galaxies appear to have ^ 0° (or, equivalently, ~ 
180°). This could be due to either: 

(i) some feature intrinsic to the code such as the initial 
estimate of ^s; or, 

(ii) a feature intrinsic to the data. 

Such a biased distribution does not arise when fitting mock 
images constructed either internally by Galactica or ex- 
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Figure 7. Histograms of parameters recovered by Galactica from our sample of SDSS galaxies, a) Characteristic sizes of the spheroidal 
and disk components, b) Ellipticity of the spheroidal components, c) Position angles of disks and spheroids. The dashed histogram shows 
the position angles of well resolved spheroids (those with effective radii greater than three times the PSF FWHM). d) Cosine of the 
inclination angle of the disks. Where appropriate, spheroids are represented in red and disks in blue. 



ternally by IRAF (see Appendix IA3.2I for details of these 
tests), suggesting that (i) is not the correct explanation. 

Explanation (i) can, in fact, be ruled out by rotating 
the images by some angle prior to fitting. If the problem 
were intrinsic to the code, we would expect to see no change 
in the distribution of ^s- In fact, when we rotate the images 
by 90°, we find that the distribution of is shifted by 90° 
(see Fig. [9]), indicating that it is some feature of the images 
themselves that is causing this problem. The same is true 
if we instead rotate galaxies by 45°. (Note that, in the case 
of a 45° rotation we crop to the largest square which fits 
within the rotated image. As a result, there are fewer pixels 
available to fit and therefore larger errors in the fit param- 



eters.) Point (ii) is a plausible explanation since the data 
were taken in drift-scan mode along the easterly direction 
which corresponds to ^ = 0°. This can lead to small asym- 
metries in the actual PSF (Berstein & Jarvis 2002). Since 
we are using a circularly symmetric PSF in our photometric 
model, Galactica may try to fit slightly elliptical bulges 
with ~ 0° to match the actual PSF shape. Note that, as 
expected, for well-resolved spheroids, the distribution of 
is close to uniform (dashed histogram in the lower-left panel 
of Fig.EI. 

For our purposes, the crucial issue is whether the bias in 
affects the derived S/T. To quantify the effects of this bias 
on the recovered S/T ratio, we re- fit a sample of our images 
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Figure 9. A comparison of the bulge position angles recovered for 
a sample of ~ 100 SDSS galaxies fit before (x-axis) and after (y- 
axis) rotation by ^rotated = 45° and 90°. When rotating images 
by 45° we crop them to the largest square which fits entirely 
within the rotated image. As a result there are fewer pixels to fit 
and therefore a larger scatter in values recovered. 

keeping equal to Od (i.e. we did not allow the spheroid 
position angle to vary freely). We find that the S/T ratios 
recovered correlate extremely well with those found with 
as a free parameter, with scatter consistent with the fitting 
uncertainties in S/T (see the left-hand panel of Fig. I1Q|) . An 
excellent correlation is also found if we rotate our images by 
45° (right hand panel of Fig. I1Q|) . The larger scatter in this 
case is caused by the reduced number of pixels available for 
fitting in our rotated images. We conclude that this bias in 
the distribution of does not affect our estimates of S/T. 
We have further found that the bias in the distribution of 
is strongest for poorly resolved, low ellipticity spheroids. For 
larger spheroids, particularly those which are quite elliptical, 
there is no apparent bias. 

In conclusion, the bias in seems to be due to some 
feature intrinsic to the data, perhaps an asymmetry in the 
PSF due to the observing method. We do not believe that 
this bias affects the recovered S/T ratios at any significant 
level since re-fitting the images with the bulge position an- 
gle locked to equal that of the disk (which is essentially 
unbiased — see Fig. [7| does not significantly alter the S/T 
ratio in the vast majority of cases. 

3.2.3 S/T vs. disk inclination 

A large number of objects in the sky which are randomly 
inclined to the line-of-sight should have a uniform distribu- 
tion of cos(i). Fig. [7] clearly shows that this is not the case 
for the inclination angles of the disk components obtained 
by decomposing our sample of SDSS galaxies. 

To test whether the apparently incorrect recovery of the 
disk inclination is an artifact of the fitting procedure, a sam- 
ple of 200 mock galaxies was generated using the Galactica 



code (see Appendix [A}. The S/T ratios were chosen at ran- 
dom in the interval [0, 1]. The remaining parameters, includ- 
ing the value of cos(i), were also chosen at random. Fig. [T3l 
demonstrates that the Galactica code generally recovers 
the cos(i) distribution for 200 model galaxies quite well. 
However, a noticeable feature is a slight deficit dX i — 90° 
and a corresponding at i ~ 75 — 80°. This excess reflects 
the fact that fits avoid the 90° limit since this would corre- 
spond to fitting an infinitely thin edge-on disk and, because 
of seeing, the disks are never infinitely thin edge-on. (The 
feature remains even if the allowed inclination range is in- 
creased from [0°,90°] to [-180°, 180°].) The S/T ratios are 
not affected by this problem, i.e. model galaxies with input 
value i ^ 90° but recovered value i ^ 85° show no bias in 
the recovered S/T ratio. 

Such biases in the distribution of cos(i) have been seen 
in other studies employing 2D galaxy decomposition tech- 
niques (see, for example, Simard et al. 2002; Tasca & White 
2006) and can occur because the fitting codes tend to fit 
disk components to radial variations in axial ratio or posi- 
tion angle in spheroids (Simard et al. 2002). Tasca & White 
(2006) used Gim2D to fit the 2D images of galaxies in the 
SDSS. They found a biased distribution of cos(i), with in- 
trinsically brighter galaxies showing the most biased distri- 
bution. Fig.fTTIreproduces Fig. 10 of Tasca & White (2006), 
with results from this work overlaid. Our results, using the 
same dataset but a different galaxy decomposition code, are 
in excellent agreement with those of Tasca & White (2006). 

Allen et al. (2006) performed 2D galaxy decompo- 
sitions, also using Gim2D, on galaxies in the Millennium 
Galaxy Catalogue and found that disk-dominated galaxies 
(S/T< 0.8) had a more uniform (although still biased) dis- 
tribution of cos(i). Fig. 1121 reproduces their results, with 
comparable results from our own work overlayed. In this 
case, we find the opposite trend: our cos(i) distribution is 
more uniform for the S/T> 0.8 sample, although the errors 
are large. We find that galaxies must have angular sizes of 
several times the seeing half-width at half-maximum in or- 
der for the inclination to be well constrained. From Fig. 1 of 
Allen et al. (2006), we would therefore conclude that a large 
fraction of their galaxies should have poorly constrained in- 
clinations. 

Assuming that the bias in cos i arises due to Galactica 
using the disk component of the photometric model to fit 
radial variations in the spheroid, it is possible to make an 
approximate correction for this bias. Such a correction was 
developed by Tasca & White (2006). In g3]we will employ 
their correction, and a similar yet more detailed correction 
to assess the impact of this bias on our results. 

3.2.4 Effects of fitting Sersic index 

We have chosen not to include the Sersic index as a free pa- 
rameter in our photometric model, instead holding it fixed 
at n = 4 (corresponding to a de Vaucouler's profile). Tasca 
& White (2006) demonstrated that n = 4 provides a good 
fit to the majority of spheroids in a magnitude limited sam- 
ple, and that there is a very good correlation between the 
values of S/T obtained using fixed n = 4 and free n fits. 
To examine this in our own data we fit a subsample of our 
galaxies allowing n to vary. In Fig. [14] we show the recovered 
S/T ratios assuming a de Vaucouler's profile (x-axis) and a 
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Figure 10. Left panel: Comparison of estimates of S/T for a sample of ~ 100 SDSS galaxies estimated from the actual image and from 
the image rotated anticlockwise by 90° (black crosses). Also shown are the results of re- fitting these images with the spheroid position 
angle forced to equal the position angle of the disk (red circles). Right panel: Comparison of estimates of S/T for the same sample 
estimated from the actual image and from the image rotated anticlockwise by 90° (black crosses) and by 45° (red circles). When rotating 
images by 45° we crop them to the largest square which fits entirely within the rotated image. As a result there are fewer pixels to fit 
and therefore a larger scatter in S/T values recovered. 



Sersic profile (y-axis). There is a good correlation between 
the results obtained using de Vaucouler's and Sersic profiles. 
This is particularly true when n > 2.5. For lower values of 
n (blue points in Fig. I14p we see some large discrepancies. 
These occur for galaxies which had a low S/T in the de Vau- 
couler's fit, but are given a high S/T when fit by a Sersic 
profile. Of course, for n ~ 1 there is no difference in our pho- 
tometric model between disks and spheroids (except for the 
fact that disks may be highly inclined to the line-of-sight 
while spheroids are limited in how elliptical they may be- 
come). It is not surprising therefore that Galactica mixes 
disk light between the two model components in such cases. 
We find that, when allowing the Sersic index to be fit as 
a free parameter the fraction of light emitted by disks (av- 
eraged over all galaxies in our sub-sample using a 1/Vmax 
weighting) decreases from 60% to 52%. This effect is very 
similar to that found by Tasca & White (2006). We consider 
this to be a lower limit on the disk light fraction since, as 
discussed above, for some galaxies a fraction of the disk light 
will have been fit by a spheroidal component with n ~ 1. 

4 LUMINOSITY AND MASS FUNCTIONS 
4.1 Introduction 

The spatial abundance of galaxies is expressed by means of 
the luminosity function (LF), defined as: 

dn(M) = (/)(M)dM, (1) 

where dn is the number density of galaxies with absolute 
magnitude in the range M, M + dM. The simplest way to 
calculate the luminosity function is using the 1 /Vmax method 



in which the number of galaxies in each individual absolute 
magnitude bin is divided by the volume of space that has 
been surveyed at that magnitude. Galaxies in any given ab- 
solute magnitude range are assumed to be uniformly dis- 
tributed in the surveyed volume which is not the case if any 
local overdensities are present. Maximum likelihood tech- 
niques circumvent this problem and provide more accurate 
estimates of the luminosity function. Here, we will employ 
the l/Vmax as well as the Stepwise Maximum Likelihood 
(SWML) non-parametric estimator (Efstathiou, Ellis & Pe- 
terson 1988) which characterizes the LF as a series of 
steps. We will also employ the STY parametric estimator 
(Sandage, Tammann & Yahil 1979), assuming a Schechter 
(1976) functional form 

,^(M) = 0.41nl0,^a0-°-"(^-"^*«"+i)exp[-10-°-*<"^-^*'], (2) 

where M* is a characteristic magnitude, a is the faint- 
end slope and 0* is the normalization. Integrating over the 
Schechter function provides an estimate of the luminosity 
density. This can also be obtained by summing up all the 
individual SWML contributions. 

Computing the spheroid and disk LFs is more compli- 
cated since there is an additional constraint to be consid- 
ered (Benson et al. 2002), namely the detectability of a 
spheroid/disk depends both on the component's apparent 
magnitude and on the corresponding S/T. This needs to be 
accounted for when constructing the luminosity function. A 
detailed discussion of the application of these methods can 
be found in Benson et al. (2002). We use exactly the same 
methods as (Benson et al. 2002) to estimate luminosity 
functions from our present dataset. 
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Figure 12. The distribution of cosi split by S/T. Histograms 
show results from Allen et al. (2006) (this figure is a reproduction 
of their Fig. 9), while points show results from this work. While 
both datasets show biased distributions of cosi, the trends with 
S/T appear to differ, with the Allen et al. (2006) dataset showing 
a more uniform distribution of cosi for galaxies with low S/T. 



Figure 11. The distribution of axial ratio h/a (equivalent to 
cos i) as a function of galaxy luminosity. Histograms show results 
from Tasca & White (2006) (this figure is a reproduction of their 
Fig. 10), while points show results from this work. Both datasets 
show that the bias in cos i occurs primarily for the most luminous 
galaxies. 

Benson et al. (2002) used a functional form for STY 
parametric fits to the spheroid and disk LPs which had a 
Schechterx exponential form. We find that the functional 
form of Benson et al. does not provide a good description of 
our larger sample of galaxies. We have been unable to find a 
suitable functional form which does provide a good descrip- 
tion and so have not performed STY fits to the spheroid and 
disk luminosity function data. 

4.2 SDSS absolute magnitudes and K+E 
corrections 

In order to estimate the luminosity function, we require 
galaxy absolute magnitudes. A galaxy at redshift with 
apparent magnitude m, has an absolute magnitude M given 
by: 

m - M = 25 + 5 logio (^l) + KE(^) (3) 

where Dl is the luminosity distance in megaparsecs and 
KE(z) is the K+E correction. 

K+E corrections for our catalogued galaxies were ob- 
tained using a code kindly provided by Carlton Baugh. 
It employs the revised isochrone stellar population syn- 
thesis models of Bruzual & Chariot (1993) to determine 
present-day galaxy luminosities. The model assumes a stel- 
lar initial mass function (IMF) and a star- format ion rate 



'0(t) (X exp(— t/r), with timescale, r. A grid of models was 
generated by varying the metallicity and r. We assume a 
Salpeter (1955) IMF and apply a simple dust extinction law. 
At every point on the grid, a table of absolute magnitudes, 
galaxy colours, K+E corrections and galaxy stellar mass- 
to- light ratio is generated. The model that best matches 
the observed g — r and r — i colours of each galaxy is then 
used to infer its present-day (z = 0) r-band absolute lumi- 
nosity, K+E correction and stellar mass-to-light ratio. The 
mass-to-light ratio is used to convert luminosities to stellar 
masses in order to estimate stellar mass functions (see ^4.5p . 
Note that the K+E corrections are based on the total (i.e. 
spheroid plus disk) colour of a galaxy. 

4.3 Luminosity function estimates 

We estimate luminosity functions using the methods de- 
scribed in detail by Benson et al. (2002) and employ both 
the SWML and 1/Vmax estimators (for the total luminosity 
function we also employ the STY estimator). We estimate 
the luminosity functions of spheroids and disks, as well as 
the total galaxy luminosity function for our sample of SDSS 
EDR galaxies. 

As noted in Appendix IB2[ our requirement that im- 
ages be reasonably well fit by Galactica (i.e. < 2) 
introduces some bias in both the apparent magnitude and 
redshift distributions of our galaxy sample. To correct this 
bias we make the assumption that the distribution of S/T 
for galaxies with > 2 is the same as that for galaxies 
of comparable apparent magnitude and redshift and with 
< 2. Such an assumption may of course not be correct, 
for example if disk-dominated galaxies are more likely to be 
poorly described by our photometric model. Nevertheless, 
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Figure 13. The input (black) and the recovered (blue) cos(i) 
distribution for 200 model galaxies created and decomposed us- 
ing the Galactica code. The figure demonstrates that the non- 
uniformity in the cos(i) is not caused by the fitting code. An 
apparent excess of galaxies with i ~ 75 — 80° can be seen. 
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this assumption represents the simplest correction that can 
be made for the bias. Therefore, for each galaxy with > 2 
we identify all well-fit galaxies with apparent magnitude dif- 
fering by less than 0.1 and redshift differing by less than 0.03 
from the true values for the poorly-fit galaxy. We then se- 
lect a galaxy from this sample at random and adopt its S/T 
ratio for our poorly-fit galaxy. 

In Appendix I A3. 2 1 we find that the value of S/T re- 
covered by Galactica (and also Gim2D) for mock images 
are biased. The median bias in S/T produced by Galac- 
tica can be approximated by a linear dependence on the 
true S/T (see Fig. IA6|) . We use this linear relation to apply 
a correction to the value of S/T recovered for each SDSS 
galaxy in order to obtain an estimate of the unbiased value. 
We use these corrected estimates of the S/T when estimat- 
ing luminosity and mass functions. 

We find that there are only small changes in the mea- 
sured luminosity functions, the most significant being a 
small enhancement in the abundance of bright spheroids. 
The luminosity density ratio quoted above varies by less 
than 0.5a after correcting for this bias. 

Our results are displayed in Fig. [15] For galaxies whose 
images are well-fit by our model, we find that the STY 
method accurately recovers the parameters of the total lu- 
minosity function; furthermore, the STY fit traces the cor- 
responding SWML points very well. The values of M* and a 
obtained from the STY fit to the total luminosity function 
agrees very well with that of Nakamura et al. (2003) (SDSS 
r-band, z = 0). While we have not been able to find a para- 
metric form which fits the spheroid and disk 2D luminosity 
functions {^{M, S/T)) we have determined the parameters 
of Schechter functions which fit the SWML data points rea- 
sonably well. These should not be considered good fits in a 
statistical sense, merely useful fitting functions. The param- 
eters of the best fit Schechter functions are given in Table [1] 

We calculate luminosity densities of disks and spheroids 
by integrating over the SWML pointfl We find the lu- 
minosity densities for spheroids and disks to be: pL = 
O.611±O.OO8xlO^/iL0 Mpc-2 and pL = L07=b0.02x lO^/iL© 
Mpc~^ respectively. 

These values are in contrast with the findings of the pre- 
vious study of Benson et al. (2002) who found the spheroid 
and disk luminosity densities to be very nearly equal. Of 
course, Benson et al. (2002) used a very small sample of 
galaxies to compute luminosity densities, finding a ratio of 
disk to spheroid luminosity density of 1.2 ±0.9. Our current 
sample gives a ratio of 1.75 ± 0.04, which is consistent with 
that of Benson et al. (2002). 



Figure 14. The S/T ratios measured for a subsample of our 
galaxies. On the x-axis we plot the value obtained assuming a de 
Vaucouler's profile for the spheroid (i.e. a Sersic profile with fixed 
n = 4) while on the y-axis we show the results of fits in which 
we allow the Sersic index n to vary. Red points show galaxies for 
which the best fit n > 4, magenta points show galaxies for which 
2.5 < n ^ 4 while blue points show galaxies with n ^ 2.5. 



4.3.1 Corrections for Systematic Effects 

Tasca & White (2006) propose a method to correct for 
the bias introduced by the non-uniform cosi distribution 

^ No correction is included for galaxies fainter than the lower 
limit shown in the figures. Using the best-fit Schechter func- 
tions listed in Table ^ we estimate that including fainter 
spheroids/disks would lead to corrections of l%/4% respectively. 
We regard these corrections as speculative since the Schechter 
function does not provide a good fit to the spheroid and disk 
luminosity functions. 
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Table 1. Best fitting Schechter function parameters for luminosity functions of total, spheroid and disk components. For the total 
luminosity function the best fit parameters are determined using the STY method. For the disk and spheroid luminosity function 
we instead fit a Schechter function to the non-parametric luminosity function determined using the SWML method — these should be 
considered useful fitting functions only, not good fits in any statistical sense. For the spheroid and disk luminosity function fits, the 
maximum deviation from the SWML data points is given in the final column. 



Component * — 5 log cx 0o/^^Mpc~^ Max. Dev. 

Total -20.62 -1.19 0.00155 N/A 

Spheroid -20.98 -1.18 0.00348 32% 

Disk -20.40 -1.39 0.00830 27% 



discussed in ^3.2.31 This method involves using only those 
galaxies with cosi < 0.5 (which Tasca & White (2006) 
consider to be true disks) to estimate the fraction of light 
emitted by disks as a function of absolute magnitude. This 
function, /disk(Mr), is then averaged over the total galaxy 
luminosity function in order to obtain an estimate of the 
fraction of light emitted by disks. Tasca & White (2006) 
include a correction for the inclination-dependent dust ex- 
tinction experienced by galaxy disks, finding that, at any 
given magnitude, /disk should be 2.56 timefl the fraction of 
light emitted by disks with cosi < 0.5. Figure [16] shows the 
function /disk for our galaxies. Filled red points are the result 
of summing the luminosities of disks over all values of cosi, 
i.e. with no attempt to correct for the biased distribution of 
cosi. Filled black points show the result after applying the 
Tasca & White (2006) correction (note that in cases where 
this correction would imply /disk > 1 we limit the value to 
unity). Averaging over the SDSS r-band luminosity function 
of Blanton et al. (2003) we find that (53 ± 3)% of the local 
luminosity density is contributed by disks. This is consistent 
with the (54 ± 2)% obtained by Tasca & White (2006). It 
should be noted that this result is robust to changes in our 
decision to include all galaxies with < 2 in our final sam- 
ple. Reducing this cut to < 1-2 for example results in a 
disk luminosity fraction of (47 =b 3)% — consistent with the 
previous result within the quoted errors. 

We can attempt to use this same approach to construct 
disk and spheroid luminosity functions corrected for the non- 
uniform cosi. To do this, we take our catalogue of galax- 
ies and identify those with cosi < 0.5 These galaxies are 
assumed to have been correctly fitted (i.e. the disk com- 
ponent of our fit corresponds to a real thin disk in these 
galaxies) and are placed into a refined catalogue. Since we 
assume that the true cosi distribution should be uniform 
we expect one galaxy with cosi > 0.5 for each galaxy with 
cosi < 0.5. Therefore, for each galaxy in our cosi < 0.5 sam- 
ple we search for a galaxy with similar spheroid and face-on 
disk absolute magnitudes but with cosi > 0.5. The most 
similar galaxy is added to our refined catalogue. At the end 
of this procedure what remains is a sample of galaxies with 
cosi > 0.5 for which there are no cosi < 0.5 counterparts. 
We assume that in these cases the disk component has been 
used to fit some feature in the spheroid. Therefore, we set 
the S/T ratio of these remaining galaxies to 1 and include 
them in our refined catalogue. 

^ This correction would be precisely 2 if there were no dust- 
extinction of the galaxy disks. 



This procedure should give a conservative lower limit 
to the disk luminosity fraction. The disk fraction obtained 
via this method is shown by the blue points in Fig. 1161 Note 
that this matches the Tasca & White (2006) method for 
bright galaxies, but falls below it at faint magnitudes. The 
reason for this is simple: the Tasca & White (2006) method 
assumes that the total disk luminosity in any bin of absolute 
magnitude is 2.56 times that of disks with cosi < 0.5 in that 
bin, even if that exceeds the total fitted disk luminosity of all 
galaxies in that bin of absolute magnitude. Thus, the Tasca 
& White (2006) method can create additional disk light in 
some bins, contrary to the assumption that the image de- 
composition code has added in extra disk light to fit details 
of the spheroid component. The open red circles in Fig. [T6l 
show the fraction of light from disks with cosi < 0.5 while 
the solid red line indicates the total disk luminosity reduced 
by a factor 2.56. Where the open red circles lie above the 
red line the Tasca & White (2006) method must create ad- 
ditional disk light. In our more detailed method, disk light 
can never be created, and so the blue points always lie below 
the red points. 

The problem just discussed illustrates the limitations 
of the Tasca & White (2006) method, and indicates that 
the reality here is significantly more complicated than the 
simple assumption adopted by Tasca & White (2006). Nev- 
ertheless, our more detailed implementation of their method 
should still give a good lower limit on the disk luminosity 
fraction. The resulting disk and spheroid luminosity func- 
tions are shown as crosses in Fig. 1151 We find a disk lumi- 
nosity fraction of (43 =b 1)%. 

In short, the Tasca & White (2006) method works pro- 
vided all objects with cosi < 0.5 are correctly fit (i.e. the 
model disk is fit to a true disk). // this assumption is correct, 
then our data imply that Galactica and Gim2D must be 
systematically failing to fit the disk components of equiva- 
lent galaxies with cos i > 0.5, assigning some of the disk light 
to a spheroid component. This could occur, for example, if 
in face on galaxies the codes use the spheroid component to 
fit a bar feature in the disk. 

To summarize, our results suggest that stars in disks 
contribute between 43 and 64% of the local luminosity den- 
sity. Tasca & White (2006) found a disk contribution of 
54it2% which is entirely consistent with this range. Fur- 
thermore, if we apply Tasca & White's correction for the 
bias in cos i precisely as they did we find a disk fraction of 
(53 ±3)%, in excellent agreement with their result. How- 
ever, as we have shown above, it is not clear that the Tasca 
& White (2006) correction is entirely valid and hence we 
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Figure 15. Luminosity functions for our sample of 7493 galaxies 
with X? < 2.0. Magnitudes are the total absolute magnitudes of 
the galaxies in the upper panel, and the absolute magnitudes of 
spheroid and disk components in the lower panels. Open symbols 
show the 1/Vmax and SWML estimates, while crosses show he 
SWML estimate after correcting for the biased distribution of 
cosi; the solid line in the top panel represents the STY fit, while 
dotted lines in the lower panels indicate the best fit Schechter 
function to the SWML data points. The top panel displays the 
total galaxy luminosity function, the middle panel the luminosity 
function of spheroids and the lower Danel, the luminosity function 



i 




Figure 16. The disk light fraction, /disk, as a function of abso- 
lute magnitude. Red points show /disk measured from our image 
decompositions with no attempt to correct for the non-uniform 
distribution of cos i. Black points show the disk fraction resulting 
from the correction described by Tasca & White (2006) (error 
bars are shown only for these points for clarity — they are simi- 
lar for other points). Blue points show the result of applying a 
more detailed correction motivated by the assumptions of Tasca 
& White (2006). Red open circles indicate the fraction of light 
in disks with cosi < 0.5, while the red line shows the filled red 
points reduced by a factor of 2.56. 



prefer to quote the range above which we feel is a very con- 
servative estimate of the disk contribution to the luminosity 
density. 

Finally, as noted in ^3.2.1[ a we suspect that Galac- 
TICA frequently uses a highly elliptical spheroid component 
to fit bar-like features in galaxy disks. If we assume that 
all spheroids at the upper limit of allowed ellipticities (i.e. 
those in the final bin in Fig. [T]) are in fact bars, and there- 
fore count their light as originating from the disk we find 
that our estimate of the disk luminosity density is increase 
by 16% while that of the spheroid luminosity density is de- 
creased by 10%. Consequently, this correction would adjust 
the disk luminosity density fraction from 43% up to 50%. 

4.4 Comparison with theoretical predictions 

In Fig.[T7]we compare our estimate of the disk and spheroid 
luminosity functions with predictions from the Baugh et al. 
(2005) and Bower et al. (2006) implementations of the gal- 
form semi-analytic model of galaxy formation. These two 
models differ in a number of important respects. For exam- 
ple, in Bower et al. (2006) feedback from the emission of 
active galactic nuclei plays a role in quenching cooling flows 
in clusters; in the Baugh et al. (2005) model, a top-heavy 
IMF is assumed for stars that form in starbursts. The two 
models, however, assume similar mechanisms for the forma- 
tion of disks and spheroids: disks form when spinning gas 
cools in a halo while spheroids form either by major merg- 
ers or by instabilities in the disks. Although both models 
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Figure 17. Luminosity functions disk and spheroid and to- 
tal light (indicated by colour; see legend for details). The sym- 
bols show our estimates for SDSS galaxies in this work and the 
lines two different implementations of the GALFORM semi-analytic 
model: Baugh et al. (2005) (dotted lines) Bower et al. (2006) 
(dashed lines). 



generally provide a reasonable description of many galaxy 
properties, they have different strengths and weaknesses. 
Neither of them has been previously applied to the study of 
galaxy morphology, although the parameter /eiiip (first in- 
troduced into semi-analytic models by Kauffmann, White & 
Guiderdoni (1993)), which controls the mass ratio at which 
a galaxy merger is deemed to destroy any pre-existing disks 
and create a spheroid, was constrained to produce a good 
match to morphological fractions as a function of absolute 
magnitude and galaxy colours. 

Fig. [17] shows that both GALFORM models reproduce 
the main trends seen in the SDSS luminosity functions. At 
faint magnitudes, the luminosity function is dominated by 
disks while at bright magnitudes disks and spheroids make 
comparable contributions. The Bower et al. (2006) model 
in particular provides a good match to the SDSS luminosity 
functions. 
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Figure 18. Stellar mass functions for galaxies as a whole and 
their disk and spheroid components (upper and lower panels) ob- 
tained using the SWML estimator. For the total mass function we 
also plot the Schechter function derived using the STY method 
(solid line), with the constraints on the parameters M* and a 
shown in the inset. For the disk and spheroid mass functions, 
the dotted lines show the Schechter function which best fits the 
SWML data points. 



4.5 Stellar mass functions 

As noted in ^4.2[ our procedure for determining K+E cor- 
rections also provides an estimate of the stellar mass of each 
galaxy. Using these stellar masses we have constructed total, 
spheroid and disk stellar mass functions using the SWML 
and STY (for total mass only) methods. For the spheroid 
and disk mass functions we also derive the Schechter func- 
tion which best fits the SWML data points. It should be 
noted that we implicitly assume that the mass-to-light ra- 
tio, T, determined for each by our K+E correction procedure 
is the same for both the disk and spheroid components. In 
reality, the recovered value of T reflects some weighted av- 
erage of the T of each component. To improve upon this sit- 
uation would require a more advanced procedure in which 



the T (and K+E correction) of disk and spheroid compo- 
nents were estimated separately using measurements of the 
disk and spheroid colours. This would require performing 
spheroid-disk decompositions in multiple wavebands. 

A quantity of interest is the ratio of stellar mass in disks 
to that in spheroids, averaged over the entire galaxy popu- 
lation. Integrating the stellar mass densities we obtain the 
average density of stars in disks and spheroids in units of the 
critical density. We find Qstars,disks = (0.486±0.004)/i"^ 10"^ 
and ^^stars, spheroids = (0.465 ± 0.006)/i~"^ 10~^ . These re- 
sults are in good agreement with those of Benson et al. 
(2002) who found atars,disks = (0.51 ± 0.08)/i"^ 10"^ and 
^stars,spheroids = (0.39 ± 0.06)/i~^ 10~^ . Wc coucludc that 
the fraction of stellar mass found in disks today is 51 ± 1%. 

If we adopt the same correction for the biased distribu- 
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Table 2. Best fitting Schechter function parameters for stellar mass functions of total, spheroid and disk components. For the total 
stellar mass function the best fit parameters are determined using the STY method. For the disk and spheroid stellar mass functions 
we instead fit a Schechter function to the non-parametric stellar mass function determined using the SWML method — these should be 
considered useful fitting functions only, not good fits in any statistical sense. For the SMBH mass function we find that a generalized 
Schechter function (see eqn. 0J provides a better fit to the data. The 7-parameter of this function is given in the final column. 



Component 


logio(M*A-^M0) 


cx 


cf)o/h^Mpc-^ 


7 


Total 


10.82 


-1.57 


0.0035±0.0002 


1 


Spheroid 


10.87 


-0.79 


0.0019 


1 


Disk 


10.64 


-0.78 


0.0035 


1 


SMBH 


7.61 


-0.65 


0.0029 


0.6 



tion of cos i as we used in ^4.3.1l we can construct stellar mass 
functions for disks and spheroids. The results are shown in 
Fig. I18[ with parameters of Schechter function fits given 
in Table (2] After applying this correction we find stellar 
mass densities of llstars,disks = (0.330 ± 0.004)/i-^ 10"^ and 
atars,spheroids = (0.622 db 0.010) /i" ^ 10"^ SO that (35 it 1)% 
of stellar mass at z = is found in disks. 

As discussed above, a difficulty in converting from light 
to stellar mass is that we expect that disks and spheroids 
should have rather different mass-to-light ratios. In the 
above, we have used a mean mass-to-light ratio, estimated 
from our K-hE method, to convert disk and spheroid light 
to disk and spheroid stellar mass. To examine the conse- 
quences of this, we perform the following simple experiment. 
We use our dataset to find the mean mass-to-light ratios as 
a function of redshift of systems identified as pure disks and 
pure spheroids (technically we identify systems which are 
at least 90% disk or 90% spheroid respectively). We then 
assume that in composite systems (i.e. galaxies with com- 
parable fractions of light in their disk and spheroid) the 
mass-to-Ught ratios of the individual components are given 
by these mean values for pure systems. We can then esti- 
mate the stellar mass of the disks and spheroids using the 
total luminosity, measured S/T and the estimated mass- 
to-light ratios for disk stars and spheroid stars separately. 
We then compute stellar mass densities using these revised 
masses. We find that this changes our results by less than 
the errors quoted above. This approach represents only an 
approximate method for finding the mass-to-light ratios of 
individual components. Nevertheless, it suggests that such 
corrections will be small. 




Figure 19. The mass function of supermassive black holes 
in galactic spheroids. Symbols show the black hole mass 
function implied by our observationally determined spheroid 
stellar mass function assuming that M^/Mq = 1.6 x 
lO^[Mgpheroid/lO^^M0]i i2 (Raring &; Rix 2004). The dotted 
line shows the Schechter function which best fits the SWML data 
points. Other lines show results from the galaxy formation model 
of Malbon et al. (2006) when using the parameters specified 
by Baugh et al. (2005; dot-dashed line) and Bower et al. (2006; 
dashed line). 



4.6 Black hole mass function 

In the last few years, it has been conclusively demonstrated 
that many galaxies posses central supermassive black holes 
and that their mass is strongly correlated with the prop- 
erties of the galaxy's spheroid such as its luminosity, stellar 
mass and velocity dispersion (Kormendy & Richstone 1995; 
Magorrian et al. 1998; Merritt k Ferrarese 2001; Marconi 
& Hunt 2003; Haring & Rix 2004). Although there is only 
direct evidence for these black holes in bright galaxies, it 
seems quite plausible that galaxies of all sizes have a central 
black hole (e.g. Malbon et al. 2006). 

From the mass function of galactic spheroids deter- 
mined in ^4.51 assuming that all spheroids contain a su- 
permassive black hole at their centre, we can estimate the 



mass function of supermassive black holes in the local Uni- 
verse. We assume that the black hole mass is given by 
M./Mq = 1.6 X lO^[Mspheroid/lO^^M0]i-i2 (Haring & Rix 

2004) and ignore any scatter in this relation since an accu- 
rate determination of the black hole mass function would 
first require a deconvolution of the (uncharacterized) error 
distribution of spheroid masses. 

The resulting black hole mass function is shown in 
Fig. 1191 Integrating this mass function gives a total black 
hole mass density in the local universe of = (2.8 ± 0.7) x 
lO^/iMoMpc"^ where we have included a scatter of 0.3 dex 
in the spheroid-SMBH mass relation (Haring & Rix 2004) 
and have included the error in the zero-point of the (Haring 
& Rix 2004) relation in our error budget. This result is 
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in agreement with previous determinations (Yu & Tremaine 
2002; Aller & Richstone 2002; Marconi et al. 2004; McLure 
& Dunlop 2004; Shankar et al. 2004) based on much smaller 
samples of galaxies. 

Applying our correction for the biased distribution of 
cos i results in an SMBH mass function shown by the crosses 
in Fig. [19] We find that a generalized Schechter function of 
the form 

*,(t)-»p[-(t)l w 

provides a better fit to this SMBH mass function. Param- 
eters of the generalized Schechter function which best fits 
the SWML data points are given in Table [2 The SMBH 
mass density after applying this correction is p, = (3.77 it 
0.97) X lO^/iMoMpc"^, which is consistent with previous 
determinations. 

For comparison with our inferred black hole mass func- 
tion, we show results from the recent model of Malbon et 
al. (2006) who incorporate a calculation of supermassive 
black hole growth into the Galform semi-analytic model of 
galaxy formation in a ACDM universe using methods similar 
to those first described by Kauffmann & Haehnelt (2000). 
The lines in Fig. 1191 show their results for two different spe- 
cific galaxy formation models. While the calculation based 
upon the parameters of Baugh et al. (2005) seems to match 
the abundance of the high-mass black holes quite well, nei- 
ther model is able to reproduce the inferred low abundance 
of less massive black holes. 

Before our results can be used to constrain such models 
strongly, it will be necessary to achieve a significantly better 
understanding of the uncertainties in the measured spheroid 
mass, and to perform the conversion from luminosity to stel- 
lar mass using a technique which accounts for the different 
stellar populations in the spheroid and disk. 

4.7 The distribution of S/T 

Finally, we examine the distributions of S/T obtained af- 
ter applying the corrections for systematic effects described 
previously in this section. Not only are these distributions 
of interest in their own right, but they can also serve as 
valuable checks for additional systematic biases in our fit- 
ting procedure. We expect that the distribution of S/T (and 
therefore the mean S/T) should be independent of disk in- 
clination and of redshift (at least for a sufficiently shallow 
sample that evolution can be ignored). In Fig. [20] we show 
the mean S/T (weighted by 1/Vmax) as a function of these 
two quantities, split by galaxy luminosity and redshift. The 
mean S/T seems to be reasonably independent of disk in- 
clination, although there is some evidence for a rise at low 
cos(i). Note that we have excluded all galaxies for which we 
believe that the disk component of the photometric model 
has been used to fit some aspects of the true spheroid. 

When we consider the mean S/T as a function of red- 
shift we see a significant increase in S/T at low-redshifts for 
the fainter samples. Our sample is not sufficiently large to be 
unaffected by large-scale structure and, in fact, shows clear 
evidence of peaks in the redshift distribution presumably 
caused by large scale structure (see Fig. IB5|) . This could 
create a redshift dependence in the mean S/T if, for ex- 
ample, a cluster of galaxies (likely to contain a substantial 



population of elliptical galaxies) is present in the sample 
at low redshift s. Larger samples, unaffected by large scale 
structure, would be needed to address this issue further. For 
now, we simply note that for our brightest cut, the mean 
S/T seems quite independent of redshift. 

Finally, we show in Fig. [21] the median S/T (and 10^^ 
and 90^^ percentiles) as a function of stellar mass. There is a 
strong trend for increasing S/T with stellar mass — the most 
massive galaxies are ellipticals. Interestingly, the median 
S/T is fairly constant at ~ 0.3 below around 3 x lO^'^/i^M©, 
after which it rises rapidly to become close to unity. This 
is similar to the 3 x 10^° found by Kauffmann et al. 
(2003) to mark the division between galaxies with young 
stellar populations, low surface densities and low concentra- 
tions and those which are older, higher density and more 
concentrated. 



5 SUMMARY AND CONCLUSIONS 

We have used a sample of ^ 9000 galaxies extracted from the 
Sloan Digital Sky Survey to estimate the spheroid and disk 
luminosity and stellar mass functions in the local universe 
using the Galactica code of Benson et al. (2002). The 2D 
model fits to the surface brightness have revealed a bias in 
the recovered disk inclination angle arising from the lack of 
strong constraints on this parameter for most galaxies. 

We find that at faint r-band luminosities, the light is 
dominated by disks whereas at bright luminosities, it is 
dominated by spheroids, with the changeover occurring at 
around the characteristic luminosity L*. Integrating the lu- 
minosity functions, we find the total r-band luminosity den- 
sities in spheroids and disks to be pL — 0.611 =b 0.008 x 
10^ /iL© Mpc-^ and pL = 1.07 ± 0.02 x 10^ /iL© Mpc"^ 
respectively. Thus, the disks contribute approximately two 
thirds of the total luminosity density. This is in contrast 
with the findings of previous studies (Schechter & Dressier 
1987; Benson et al. 2002), based upon galaxy samples over 
40 times smaller, which found the spheroid and disk lumi- 
nosities to be very nearly equal. 

Due to the fact that real galaxies do not always look 
like our idealized models we find a biased distribution of disk 
inclinations cosi. This bias has been noted before by Simard 
etal. (2002), Tasca& White (2006) and Allen et al. (2006). 
The bias found here is identical to that found by Tasca & 
White (2006). Attempting to correct for this bias leads us 
to a conservative estimate for the disk contribution to the 
local luminosity density of between 43 to 64%. Applying the 
correction suggested by Tasca & White (2006) we find a 
disk contribution of (53 it 3)% in excellent agreement with 
their result of (54 ±2)%. 

Current a priori galaxy formation models are able to 
reproduce the disk and spheroid luminosity functions rea- 
sonably well — at least to the extent of predicting the correct 
trends of abundance with luminosity. 

Using an approximate conversion of r-band light to stel- 
lar mass, we derive ^2stars,disks = (0.486 db 0.004)/i~^10~^ 
and r^stars, spheroids = (0.465 ± 0.006)/i~^ 10~^ , in excellent 
agreement with the earlier work of Benson et al. (2002). 
Correcting for the bias in the cosi distribution leads to re- 
vised values of ^^stars, disks = (0.458 ± 0.005)/i~^ 10~^ and 
^stars,spheroids = (0.622zb0.010)/i"^ 10"^. Thus, stars in disks 
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Figure 20. The mean, l/Vmax weighted S/T ratio as a function of disk inchnation (left-hand panel) and redshift (right-hand panel). 
Points show the mean S/T in each bin, while error bars show the Icr error on the mean. In the left-hand panel, galaxies are split into four 
groups by luminosity and redshift at Mr — Slog^^Q h < —20.57 and z < 0.085. The bright, low-redshift sample is shown by red points, the 
bright high-redshift sample by blue points, the faint low-redshift sample by magenta points and the faint high-redshift sample by green 
points. Galaxies for which the disk component of the model is thought to have been used to fit features of the observed spheroid have 
been excluded. In the right-hand panel galaxies are split by luminosity as indicated by the labels in the figure. 
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Figure 21. The distribution of S/T as a function of galaxy stellar 
mass. Galaxies are weighted by l/T/max- Note that many galaxies 
are assigned S/T= 1 during our process of correcting for cases 
where the disk in the photometric model has been fit to a true 
spheroid component. Red points with error bars show the median 
S/T as a function of stellar mass together with the 10^^ and 90*^ 
percentiles of the distribution. 

contribute between 35 and 51% of the local stellar mass den- 
sity. Bell et al. (2003) claim that between 50 and 75% of 
the stellar mass density comes from late-type galaxies. This 
range is in excellent agreement with our results. It should be 



noted that the morphological selection chosen by Bell et al. 
(2003) would essentially select all galaxies with a spheroid- 
fraction of 0.5 or greater (and a large number of galaxies 
with smaller spheroid fractions) according to our calcula- 
tions (see Fig. IC2|) . 

From the inferred spheroid mass function and the ob- 
served relation between central supermassive black hole 
mass and spheroid stellar mass, and assuming that all 
spheroids harbour a central black hole, we infer the su- 
permassive black hole mass function. The associated black 
hole mass density in the local universe is p« = (3.77 ± 
0.97) X IO^/iMqMpc""^, consistent with previous estimates. 
Improvements in the characterization of the errors in the 
spheroid mass function and in the stellar population mod- 
elling will enable a better estimate of the spheroid (and black 
hole) mass function. 

We conclude that the local Universe contains around 
roughly comparable amounts of stars, by mass, in disks and 
in spheroids. This fundamental ratio is the outcome of the 
physical processes at play in the formation of the galaxy 
population. 
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APPENDIX A: METHODS FOR 
QUANTITATIVE GALAXY MORPHOLOGY 

In this Appendix, two independent methods for spheroid-to- 
disk decompositions of galaxies are described and compared. 
A number of tests are performed to reveal and estimate any 
potential biases in the decomposition codes. 

Al Introduction to spheroid-to-disk 
decomposition 

To first approximation, the surface brightness of a galaxy 
can be expressed as the sum of a highly concentrated central 
component, the spheroid, and an extended disk. Empirically, 
the surface brightness profiles of spheroids and disks are well 
represented by the following functions, 

/s = /eexp(-7.67[(r/re)'/^ - 1]), (Al) 

an r^/^-law, for the spheroid (de Vaucouleurs 1961), where 
Te is the half-light radius and le is the surface brightness at 
re, and, 

/d = /oexp(-r/rd), (A2) 

an exponential- law for a face-on disk, where rd is the expo- 
nential disk scale-length and Iq is the central intensity. 

Equations lAll and IA2I can be used to construct model 
images of the galaxy. Comparison with the surface bright- 
ness distribution of each galaxy, including the effects of in- 
clination, enables the fitting parameters to be determined. 

Andredakis et al. (1995) used this technique to fit the 
spheroid components of a sample of morphologically selected 
galaxies with types ranging from SO to Sbc. They assumed 
a more general type of profile, 

/s = /eexp(-6n[(r/re)'/"-l], (A3) 



first proposed by Sersic (1968), where n is often referred 
to as the Sersic index and determines the 'peakiness' of the 
profile and bn is a constant dependent on the value of n. 
Andredakis et al. found that the value of n varied system- 
atically from 1 for late-type spheroids to 6 for early-type 
spheroids, de Jong (1996) also suggested that the spheroids 
of field spirals are better fit using a pure exponential (i.e. 
n = 1) profile. 



A2 Methods for 2-dimensional spheroid-to-disk 
decomposition 

A2A Fitting parameters revisited 

On a Cartesian grid {x,y) the effective r in eqns. ()A1|) and 
(|A2)) become 

r'^{x,y) = — [xcos(^s) - 2/sin(^s)]^ + es[xsin(^s) 

+ycos{e,)f (A4) 

(A5) 

and 

r^{x,y) = [xcos(6'd) -ysm(Od)f ^ ^-^ [x sin(6'd) 

cos(zj^ 

+2/cos(^d)]' (A6) 

respectively. Here and Od are the spheroid and disk posi- 
tion angles, where a position angle is defined as the angle of 
orientation of the galaxy's main axis with respect to some 
coordinate system, and Cs is the spheroid ellipticity used to 
describe the deviation from circularity of the spheroid com- 
ponent 

In terms of a r^^^ spheroid and an exponential disk, 
and including the sky background, 2D decomposition usu- 
ally requires a total of 13 free parameters: 

• total flux in the galaxy; 

• S/T: ratio of the amount of light in the spheroid to the 
total amount of light; 

• re', effective radius of the spheroid; 

• Cs'. spheroid ellipticity; 

• spheroid position angle; 

• rd' scale length of the disk; 

• i: inclination angle of the disk; 

• ^d: disk position angle; 

• Xc,yc' subpixel offset of the galaxy centre; 

• residual sky background level; 

• FWHM of the point spread function (PSF); 

• n: Sersic index. 

In order to make the decomposition procedure as accu- 
rate and fast as possible, the following points must be taken 
into account: 

• The PSF smooths the galaxy image and so to achieve 
accurate fits, it must be modelled accurately and included 
in the mock images. 

• The fit is carried out using small "thumbnail" regions 
around each galaxy. The size of the thumbnail must be small 
enough to enable a fast fit to the image, but large enough 
to include all regions of the galaxy with significant signal- 
to-noise. 



Luminosity and stellar mass functions of disks and spheroids in the SDSS and the SMBH mass function 



• The mean sky background level should be ~ since the 
decomposition codes are designed to work with no (or very 
little) background; any excess sky light can be mistaken for 
galaxy light and lead to incorrect parameter estimation. 

We now explore the similarities and differences of two 
independent multi-dimensional fitting codes, the Gim2D 
code of Simard et al. (2002) and the Galactica code of 
Benson et al. (2002). 

A2.2 Galaxy image 2D decomposition (Gim2D) 

Gim2D is a publicly available code written by Simard et 
al. (2002) and widely used for automated spheroid-to-disk 
decompositions of galaxy light profiles (Balogh et al. 2002; 
Nelson et al. 2002). This code was purposely written for 
imaging with the Hubble Space Telescope Wide Field and 
Planetary Camera which has a very well modelled PSF 
(Krist 1995). The code can also be used for ground-based 
imaging data but, in this case, special attention must be 
given to the much larger and less well defined PSF. 

Object Detection 

To extract a postage-stamp image around each galaxy, 
Gim2D relies upon SExtractor (Bertin & Arnouts 1996). 
S Extractor determines the galaxy centroid position and 
the area at the faintest detected isophote to be obtained and 
also measures the mean level of the sky background for each 
galaxy (a 3cr threshold is usually sufficient to discriminate 
between the object and the background). Gim2D extracts 
a postage-stamp of size equal to a multiple of a galaxy 
isophotal area. A value of 15xiso_area was found to be 
optimum. The sky-background is not recommended to 
be treated as a free fitting parameter in Gim2D because 
the underlying sky is not well known and can potentially 
bias the output (Simard et al. 2002). However, before the 
decomposition procedure is initiated, Gim2D uses the pixels 
flagged by SExtractor as belonging to the background 
(flag value 0) to recompute the background value, therefore 
ensuring that the mean sky level is close to zero. All the 
background pixels and also pixels flagged as 'bad' (flag 
value -2) by SExtractor are subsequently excluded from 
the fitting altogether. 

Point Spread Function 

During the minimization in Gim2D, the PSF is kept fixed. 
For ground-based imaging, the PSF is obtained from a 
bright unsaturated stellar image; for HST data, an analytic 
PSF modelled using the Tiny Tim software is used (Krist 
1995). 

Minimization Technique: Metropolis Algorithm 
Gim2D allows up to 12 parameters to be fit|3 and uses the 
Metropolis algorithm (Metropolis et al. 1953) to search 
for the minimum in this multi-dimensional parameter 
space. Before starting the Metropolis algorithm, Gim2D 
works in the Initial Condition Filter (ICF) mode, i.e. it 
creates a user-specified number of models between the 
limits of the user-specified multi-dimensional parameter 

^ Sersic index is held fixed at 4. 



space. The ICF computes the model likelihoods and sets the 
sampling origin to the parameters of the best model, making 
it a subvolume to be exploited by the Metropolis Algorithm. 

Gim2D Outputs 

After finding the model that corresponds to the highest like- 
lihood, Gim2D produces a residual (object - model) map 
and calculates the value of the corresponding xi.li xi 1 
and the residual map has no remaining galaxy structure, the 
best-fit model is accepted. 

A2.3 Galactica 
Introduction 

The 2-dimensional decomposition code described here 
is based on a technique proposed by Wadadekar et al. 
(1999). Galactica was developed by Benson et al. (2002) 
and assumes the standard empirical formalisms for the 
2-dimensional surface brightnesses of a galaxy spheroid and 
disk components respectively f eqns. IATI IA2|) . 

Object detection 

To locate and extract a postage-stamp image around every 
galaxy, GALACTICA employs the same method as Gim2D. 
SExtractor is also used to measure the mean level of 
the sky background for each galaxy (a 3cr threshold is 
usually sufficient to discriminate between the object and 
the background) and this value is subtracted from the 
corresponding galaxy image. To mask any overlapping 
objects within the extracted postage-stamp Galactica 
relies upon an in-built masking algorithm which finds 
any objects that contaminate the galaxy of interest and 
masks them out. The galaxy itself is also detected by the 
algorithm using a 5cr threshold above the sky background. 
Pixels which have not been fiagged as belonging to any of 
the detected objects are used in the sky background fitting. 

Point spread function: Moffat profile 

To correct for the effect of seeing, the Galactica code gen- 
erates a Moffat profile star image (Moffat 1969) of a given 
full-width half-maximum (FWHM) expressed in terms of 
cpsF = PSFfwhm/2.35. This analytic profile is defined by 

PSF(r) = const/[l + (r/a) Y (A7) 

and is thought to represent the overall PSF shape better 
than a pure Gaussian which only approximates the core 
regions. Here a represents the width of the PSF and is re- 
lated to the FWHM= 2aV2i//3 _ i (Trujillo et al. 2001). (3 
governs how "peaky" the PSF profile is (the larger [3 is, the 
more Gaussian- like the profile becomes). A value of /3 = 4.5 
is used throughout. The a parameter can be fine-tuned to 
a particular dataset using the average FWHM for the data. 
Galactica lets apsp be a free fitting parameter to allow 
for any small changes in the PSF between the position of 
the stars in the image and the galaxy position. 

Minimization technique: Powell's method 
The code requires explicit initialization of the fitting param- 
eters. The initial value of the S/T ratio is always fixed at 
0.5, although starting with a value of S/T randomly dis- 
tributed in the range to 1 has no effect on the recovered 
S/T distribution. The starting position angles of the disk 
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and spheroid components, their characteristic radii and the 
disk inchnation angle are calculated directly from the image. 
The ranges over which parameters are allowed to vary dur- 
ing the fitting are specified and fitting outside these limits 
is not possible. 

is minimized in a 12-parameter space; the Sersic 
index is set to n = 4 and the FWHM of the PSF and 
the residual sky background level are additional fitting 
parameters not included in Gim2D. Note that in fitting 
the total flux in the galaxy Galactica we find a typical 
variation of only 5% around the value estimated from 
SDSS photometry. The minimization routine is also rather 
different from the one used by Gim2D. In Gim2D every 
parameter is varied at each step according to the 'tempera- 
ture' of the fit. In Galactica one parameter is minimized 
at a time, i.e. all but one parameters are 'frozen' until a 
minimum for this parameter is found and the process is 
repeated for the entire set of parameters until the global 
minimum is found — the essence of Powell's method (see 
Press et al. 1992 for further details). This method is good 
for finding a global extremum but is typically slower than 
the Metropolis algorithm employed by Gim2D. 

Galactica outputs and error estimation 
After convergence is achieved, the best-fit parameters are 
output along with the best-fit model image and the residual 
map obtained by subtracting the model galaxy from the real 
image. The value of per degree of freedom, is then 
calculated. Errors on the fitted parameters are obtained 
using a Monte Carlo method: we create 30 realizations of 
the best fitting model for each galaxy by adding random 
noise, drawn from a Poisson distribution, to the model 
image. The distribution of the best-fitting parameters of the 
model realizations is then used to estimate the uncertainty 
in the fit. This method allows the uncertainties in the image 
parameters to be obtained without any assumptions about 
their distribution. 



A3 Gim2D vs. Galactica comparison 
A3.1 Introduction 

The previous section described two independent codes for 
estimating basic galaxy structural parameters. Although 
these codes assume the same analytic surface brightness 
profiles to fit the spheroid and disk components, the 
differences in the number of fitting parameters and in the 
minimization techniques are suflficiently large to make a 
comparison interesting and important. 

We first note the following points relevant to the com- 
parison: 

• both codes assume a fixed value of n = 4; 

• the sky background is always kept fixed by Gim2D al- 
though the code is allowed to recompute and correct the 
background level before the minimization procedure starts; 

• Galactica always treats the sky background as a free 
parameter; 

• the ellipticity is defined differently: Gim2D fits e = l — 
hi a while Galactica fits a/6; 




Figure Al. Correlation between the input S/T ratios for 100 
model galaxies and the best- fit S/T ratios recovered using Gim2D. 
The scatter around the mean is <Trms = 0.10. Note a small sys- 
tematic underestimate of S/T, accompanied by increased scat- 
ter, for large input values. The recovered characteristic radii of 
these galaxies are the largest and the sky background recalculated 
by Gim2D for these galaxies is high. The resulting confusion be- 
tween the background and an extended surface brightness profile 
accounts for this effect. 
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Figure A2. Correlation between the input S/T ratios for 100 
model galaxies and the best-fit S/T ratios recovered using Galac- 
tica. The scatter around the mean is <Jrms = 0.11. 
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Table Al. Parameter ranges used for constructing mock galaxies. 



Parameter Low limit High limit 



S/T 


0.0 


1.0 


re,d (pixels) 


1 


12 


es 


0.0 


0.8 


i (degrees) 


0.0 


90.0 


^s,d (degrees) 


0.0 


180.0 


FWHM (") 


1.4 


1.4 



• the seeing is fixed in Gim2D but in Galactica it is 
allowed to fluctuate between ±5% of the specified ctpsf; 

• the position angles in Gim2D are defined with respect 
to the 2/-axis of a Cartesian system while in Galactica they 
are defined clockwise from the x-axis. (The position angles 
of the spheroid and disk are allowed to vary in both codes; a 
large difference between these can be a signature of a barred 
structure; Simard et al. 2002). 

To quantify the performance of the codes, a series of 
tests were conducted as we now describe. 

A3. 2 Tests using model galaxies 

A useful in-built feature of both Galactica and Gim2D 
is the ability to create model galaxies. The initial tests 
and code comparisons described below were carried out on 
model galaxies generated "internally" by Galactica. 

Model galaxies were constructed adopting parameter 
values chosen at random between realistic limits (Table [Al]) 
and matching the total counts measured in a typical real 
galaxy. Poisson noise was added to the model galaxy af- 
ter its image was convolved with an analytic Moffat PSF 
corresponding to a typical value of the seeing. This PSF is 
subsequently used as the Gim2D PSF. Model galaxies were 
analyzed with both codes using exactly the same procedures 
as for real galaxies. 

Comparison of the known input S/T values and the val- 
ues output by Gim2D and Galactica for 100 model galax- 
ies are shown in Figs. |AT] and [A2] respectively. In both cases 
the codes recover the input S/T very well. The scatter in 
the recovered S/T ratios IS (Trms 

^ 0.10. The remaining pa- 
rameter recoveries are shown in Figs. IA3l and |A4I 

The results shown in Figs. lAll to lA4l demonstrate that 
both Gim2D and Galactica produce, on the whole, reliable 
spheroid/disk decompositions for a set of artificial galaxies 
generated according to the model assumed by the code. A 
more stringent test of the codes is to apply them to model 
galaxies generated independently of the codes themselves. 
Model galaxies were therefore externally created using the 
IRAF task MKOBJ. 

The model parameters were taken from Table lAll a 
model galaxy is created for given values of the size, orienta- 
tion and ellipticity (in this case defined asb/a) and the image 
is convolved with a specified seeing. A useful feature of this 
approach is that a real science frame can be fed into MKOBJ 
and the model galaxy added to a blank patch of the sky on 
this science frame. By matching, on average, the counts in a 
real galaxy, this procedure ensures that the artificial image 
has similar noise characteristics to the real data. 



To generate the model galaxies several SDSS r'-band 
frames were extracted, each typically containing ^ 5 
SDSS catalogued galaxies. Each frame was taken from a 
different patch of the sky. Counts associated with the SDSS 
galaxies were measured using the S Extractor f lux_best 
parameter. In order to test the decomposition algorithm 
over a realistic range of galaxy properties, this procedure 
was applied to galaxies spanning a range of apparent mag- 
nitude and apparent shape and size. The model galaxies 
were inserted across the blank regions of the sky in the 
original frames. The post age- stamps for these galaxies were 
extracted using SExtractor and the decomposition codes 
run treating the extracted model galaxies just as the real 
ones. 

The results of the Gim2D decompositions of the model 
galaxies are shown in Fig IA5I The agreement between the 
input and output S/T ratios for the pure exponential disks 
(S/T= 0) is excellent. However, the recovered S/T ratio for 
S/T= 0.5 is biased by AS/T= 0.1 and for S/T= 1.0 it is 
biased by AS/T= 0.2. The tendency is always to underesti- 
mate the amount of spheroid or, equivalently, overestimate 
the amount of disk. Gim2D can be fine-tuned to recover the 
input S/T ratios with AS/T~ 0.1 across the full S/T range. 
For this, Gim2D requires that the size of the zone around 
the lowest SExtractor isophote used in the re-calibration 
of the sky background should be set to ^ 30 pixels (the de- 
fault value is 10 pixels). This ensures that any faint galaxy 
fiux does not contribute to the re-calibrated background fiux 
thus minimizing any bias in S/T. 

The results of the Galactica decompositions of the 
model galaxies are shown in Fig. IA6I For the pure expo- 
nential disks the recovered S/T ratios are, again, very good. 
As before, for larger S/T, this ratio is underestimated and 
peaks at S/T= 1.0, implying that many pure r^^^ galaxies 
have acquired a fictitious disk component. There appears 
to be no correlation between the size of the underestima- 
tion and the magnitude or scale radius of the input galaxy 
but a weak correlation with the minor/major axis ratios: 
the S/T deviation is largest for the most elliptical profiles. 
The most prominent correlation, however, is between the 
recovered S/T ratio and the sky background. As discussed 
earlier, the Galactica code allows the sky background to 
fiuctuate a little to allow for uncertainties in the estimated 
background. The fact that the deviation between the input 
and the output S/T ratios is largest when the 'fitted' back- 
ground is smallest implies that an extraneous disk compo- 
nent is found where, in fact, the extra counts are due to the 
sky background. Since there is no sharp cut-off in either the 
spheroid or the disk, at the galaxy edges the galaxy surface 
brightness profile and the sky background are indistinguish- 
able. 

Gim2D performs marginally better than Galactica in 
the recovery of the S/T ratios. However, both codes show 
similar biases in the decompositions. This is despite the fact 
that we allow the background and PSF to be fit by Galac- 
tica but not by Gim2D indicating that this choice does not 
significantly bias our results. The relative performance of 
these two codes on a set of real galaxies is compared next. 
This allows a more realistic comparison of the codes but, of 
course, there is no a priori correct answer. 
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Figure A3. Correlations between input and recovered values of the characteristic radius, position angle, ellipticity and inclination for 
100 model galaxies created using the Galactica code and decomposed using Gim2D. a) Disk and spheroid radii, b) Disk and spheroid 
position angles, c) Spheroid ellipticity. d) Disk inclination. There is a saturation at i = 85° which is the upper limit that Gim2D allows 
for the disk inclination. All other parameters correlate well although significant scatter is seen. 



A3. 3 Tests using real galaxies 

To ensure a uniform sampling of the [S/T, apparent 
magnitude] space, the comparison was carried out using 
a subsample of SDSS galaxies selected in bins of 0.5 in 
apparent magnitude and 0.2 in S/T ratio (as determined 
by Galactica). Unsaturated stellar images with high 
S/N were extracted from the SDSS galaxy frames and 
used in the Gim2D PSF deconvolution. The Galactica 
Moffat PSF was fine-tuned to fit the SDSS data well. 
Fig. IA7I demonstrates a significant correlation (Spearman 
rank correlation coefficient of 0.74) between the S/T 
ratios for 350 SDSS galaxies inferred using Gim2D and 
Galactica. There are no systematic differences between 



the results of the two codes. The correlations between other 
parameters are displayed in Fig. IA8I 

APPENDIX B: GALACTICA METHODOLOGY 

Bl Code speed limitations and galaxy binning 

The SDSS galaxy sample contains less than 200 galaxies 
whose post age- stamp size exceeds 91 pixels on a side. To 
reduce the processing time, these large post age- stamps were 
binned 2x2. This has the added advantage that large nearby 
bright galaxies are sampled at a resolution comparable to 
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Figure A4. Correlations between input and recovered values of the characteristic radius, position angle, ellipticity and inclination for 
100 model galaxies created and decomposed using the Galactica code, a) Disk and spheroid radii. Note that the reconstruction seems 
to hit the upper limit on the characteristic radii when the input radii are very small or very large, b) Disk and spheroid position angles, c) 
Spheroid ellipticity. Note that a number of very elliptical spheroids are found. These are mostly galaxies which have very small spheroids, 
d) Disk inclination. Note an apparent saturation at i ~ 85°. Even though the Galactica code allows the disk to be fully inclined 
(i = 90°), the reconstruction avoids this upper limit. 



that of more distant objects. To make sure that the bin- 
ning procedure does not bias the recovery of the galaxy S/T 
ratios, we carried out a series of tests. 

Several model galaxies were created using our standard 
procedure (see Appendix IA3.2|) and binned using the iraf 
task BLKAVG. GALACTICA was then used to perform the fit- 
ting ensuring that the pixel for the binned image is set to 2 x 
the normal pixel size (2 x 0.396^^) and that the noise proper- 
ties in this 'super-pixel' are changed accordingly. Figs. IBll 
and IB 2 1 show fits to model galaxies consisting either of a pure 
exponential disk or a pure r^^^ spheroid after the original 



model images were binned by 2 x 2. In the case of the pure 
exponential disk, the fit to the model is perfect. The fit to 
the pure r^^^ galaxy, however, shows a similar bias to that 
seen for the unbinned data (Appendix [A|, as indicated by 
the low value of the recovered S/T= 0.8. This shows that 
the binning in itself is not responsible for the observed S/T 
bias. The 2x2 binning was therefore applied to all the SDSS 
galaxies whose post age- stamps are greater than 91 x 91 pix- 
els. 

The binning works very well if the binned galaxy does 
not exhibit much internal structure (as in the model galax- 
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Figure A8. Correlation between various parameters inferred using Gim2D and Galactica for a sample of 350 SDSS galaxies. No 
systematic differences in the recovered parameters are apparent although the scatter can be quite large. 



ies). However, for a galaxy which exhibits significant internal 
structure, a fit with xl > 2.0 is more typical. Whether de- 
composing such galaxies even without binning would lead 
to a good fit is unclear as demonstrated for 2 SDSS galaxies 
in Fig. IB3l and lB4l The top images in both figures show the 
postage stamp and the residual map for the unbinned galaxy 
which has size 101 x 101 pixels. (To speed up the calculation, 
the postage stamp was trimmed by 5 pixels on either side.) 
The bottom panels show the corresponding images for the 
binned versions of the same galaxy. 

The galaxy in Fig. IB3I exhibits much more internal 
structure than the galaxy in Fig. IB4[ as is clear from both 
the value of xt ^ind the residual image. This supports the 
conclusion that galaxies which exhibit internal structure are 
poorly fit irrespective of whether they are binned or not. 



The recovered S/T ratios for the unbinned and binned data, 
although different, are consistent with the typical errors in 
the S/T ratio. We conclude that the S/T distribution of 
the final SDSS galaxy sample is not biased by binning these 
large, bright nearby objects (most of which contribute to the 
faint-end of the luminosity function - see Section 2]). 

B2 SDSS data and the goodness-of-fit 

The selected SDSS sample of 8839 galaxies is too large for 
each of the residual images to be inspected by eye to ensure 
a satisfactory decomposition as suggested by the xl < 2.0. 
However, a randomly selected sample of residuals was ex- 
amined by eye to ensure that they were indeed predomi- 
nately noise dominated. The xt < 2.0 criterion was therefore 
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Figure A5. The difference between input and recovered values 
of S/T using Gim2D for a set of 250 model galaxies. The mean 
offset in the recovered values are AS/T= 0.02 (for S/T= 0.0), 
AS/T= 0.13 (for S/T= 0.5) and AS/T= 0.20 (for S/T= 1.0) 
Model galaxies span a range of apparent magnitudes, sizes and 
orientations. 




Figure A6. The difference between input and recovered values of 
S/T using Galactica for a set of 250 model galaxies. The mean 
offset in the recovered values are AS/T= 0.05 (for S/T= 0.0), 
AS/T= 0.11 (for S/T= 0.5) and AS/T= 0.24 (for S/T= 1.0) 
Model galaxies span a range of apparent magnitudes, sizes and 
orientations. The median AS/T is well fit by the relation AS/T= 

0.02 -0.26S/Ttrue. 




Figure A7. The correlation between the S/T ratios of a sam- 
ple of 350 SDSS galaxies inferred using Gim2D and Galac- 
tic A. The Spearman rank correlation coefficient of 0.74 indicates 
a significant correlation. Note that almost no pure disks (i.e. 
S/T= 0.0) are detected by Galactica but a few are detected 
by Gim2D. Most of these galaxies have spheroid characteristic 
radii of less than 2 pixels (as does the single Galactica detec- 
tion at S/T= 1.0). Gim2D finds larger characteristic radii for 
these galaxies. The overall scatter is crrms = 0.19, but there is a 
marked increase in the scatter for larger S/T ratios. 



adopted to define a 'well fit' dataset of 7493 galaxies. To test 
for any selection biases were introduced by the rejection of 
galaxies with xl > 2.0, we compare the distributions of some 
basic properties of these galaxies and of the well-fit subset 
in Fig. IB5I There are small but noticeable biases introduced 
in the distributions of apparent magnitudes and redshifts. 
These biases are taken into account when estimating galaxy 
luminosity functions (see ^4.3|) The figure also shows a deficit 
of objects with S/T > 0.7. This is most likely due to the bias 
in the Galactica code discussed in Appendix [A] The sig- 
nificance of this bias and its influence on the final results is 
discussed in the main body of the paper. 

Finally, in Fig. IB6l we plot the distribution of i?50 (the 
radius enclosing 50% of the Petrosian flux) for galaxies meet- 
ing the selection criteria of this work and that of Tasca & 
White (2006). Our galaxies are typically 40% smaller than 
those of Tasca & White (2006). 

B3 S/T error estimates 

Benson et al. (2002) developed a Monte Carlo approach to 
estimate the errors on the fitted parameters in Galactica. 
This method is very time consuming, requiring several CPU 
days for a typical galaxy. A full Monte Carlo analysis is 
therefore impractical for a large dataset such as the one in 
this paper. Instead, to obtain representative error estimates, 
we split the sample into bins of apparent magnitude of width 
0.5, and selected from each of these bins 5 galaxies from each 
of three further bins in S/T (0.0 <S/T< 0.3, 0.3 <S/T< 
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Figure B5. Distributions of (a) apparent magnitude, (b) absolute magnitude, (c) redshift and (d) recovered S/T ratios. Results for the 
total sample are shown by the black histograms and for the sample with < 2.0 by the blue histograms. The two distributions are 
similar, indicating that excluding poorly fit galaxies does not introduce any obvious biases in the sample. 



0.6, 0.6 <S/T< 1.0). The full Monte Carlo analysis was 
performed on the selected subsample and the medians of 
the derived errors taken to be representative for galaxies in 
each [rmag, S/T] bin. 



APPENDIX C: SPHEROID-TO-DISK RATIOS 
AND GALAXY MORPHOLOGIES 

CI Morphological classification using colour 

It has long been known that galaxy colour is a useful in- 
dicator of whether a galaxy is elliptical (old, red) or spiral 
(young, blue) (de Vaucouleurs 1961) since the dominant 
stellar populations are reflected in the galaxy colours. Inves- 



tigating the colour-magnitude and colour-colour diagrams, 
Strateva et al. (2001) have shown that the {u — r) colour 
distribution of SDSS galaxies has two maxima which are sep- 
arated by a well-defined minimum dX {u — r) — 2.2 and that 
98% of galaxies spectroscopically classified as 'early' types 
have {u — r) > 2.2 whilst 73% of spectroscopically classified 
'late' types have (u — r) < 2.2. Strateva et al. (2001) have 
also shown that this separator also applies for a subsam- 
ple of visually classified morphological types where 80% of 
galaxies visually classified as E, SO or S a have colours redder 
than {u — r) = 2.2 and 66% of galaxies visually classified as 
Sb, Sc and Irr have colours bluer than {u — r) = 2.2. The 
{u — r) separator has already been used to study morpholog- 
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Figure Bl. Fit to a pure exponential model galaxy after the orig- 
inal 101 X 101 model image was binned 2x2. The top row shows 
the real (left) and model (right) images. The middle row shows 
the disk (left) and spheroidal (right) components. The bottom 
row shows the real (left) and residual (right) images. The recov- 
ered S/T= 0, corresponding to a pure exponential. The good fit 
is evident from both the ~ 1 and the noise-dominated residual 
image. 




S/T= 0.79 

Figure B2. Fit to a pure r-*^/^ model galaxy after the original 
101 X 101 model image was binned 2x2. The top row shows the 
real (left) and model (right) images. The middle row shows the 
disk (left) and spheroidal (right) components. The bottom row 
shows the real (left) and residual (right) images. The recovered 
S/T = 0.8, showing once more that the Galactica code returns a 
biased estimate of S/T. The good fit is inferred from both ~ 1 
and the noise-dominated residual image. 



ical properties of galaxies in the SDSS sample as a function 
of environment (Goto et al. 2002; Balogh et al. 2004). 

Our derived S/T ratios are plotted against u— r colour 
in Fig. ICll The B/T distributions of red galaxies and blue 
galaxies are significantly displaced relative to each other: 
most galaxies classified by colour as early types {u — r> 2.2) 
have S/T > 0.4 whereas most galaxies classified by colour 
as late types {u — r < 2.2) have S/T < 0.2. However, the 
red galaxies in particular span a large range in S/T ratio. 
This suggests that blue star-forming galaxies are predomi- 
nantly disk-dominated but that disk-dominated galaxies in- 
clude both star-forming (blue) and passive (red) galaxie^ 

C2 Morphological classification using 
concentration index 

Galaxies can also be classified according to how 'peaky' their 
light distribution is by using the concentration index (Abra- 
ham et al. 1994). The surface brightness distribution of 

^ The red disk population may also include galaxies with heavily 
obscured star formation. 



ellipticals and SOs is considerably more centrally concen- 
trated than that of spirals and irregulars. Shimasaku et al. 
(2001) defined the (inverse) concentration index for SDSS 
galaxies as the ratio of the half- to the 90% light radii and 
define an optimum division between late and early types to 
be at C = 0.33 (with 15 — 20% contamination from oppo- 
site types) . This separator has also been used to investigate 
the morphological properties of SDSS galaxies (Goto et al. 
2002; Nakamura et al. 2003). 

Fig. IC2I shows that there is a good correlation be- 
tween the S/T ratios and the (inverse) concentration index, 
C, for galaxies with small S/T < 0.3. Although spheroid- 
dominated galaxies tend to be more centrally concentrated 
than disk-dominated galaxies, there is little correlation of 
concentration with S/T for S/T > 0.3. There is, however, 
considerable scatter in C for a given S/T. 

C3 Morphological classification: S/T vs 
eye- morphology 

Shimasaku et al. (2001) used a sample of 456 bright SDSS 
galaxies (g' < 16.0) visually classified into seven morphologi- 
cal types (Hubble types E, SO, Sa, Sb, Sc, Sdm and Im) to in- 
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Figure B3. Fits to a galaxy that exhibits internal structure. 
The top images show the unbinned galaxy post age- stamp (left) 
and the corresponding residual image (right). The bottom images 
show the galaxy and the residual after the galaxy is binned 2x2. 
In both cases the is poor (x1 > 2.0) and the residuals are not 
noise-dominated. This supports the conclusion that galaxies with 
significant structure give poor fits irrespective of whether they 
are binned or not. 
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Figure B4. Fits to a galaxy that does not exhibit internal struc- 
ture. The top images show the unbinned galaxy postage-stamp 
(left) and the corresponding residual image (right). The bottom 
images show the galaxy and the residual after the galaxy is binned 
2 X 2. In both cases, the xi is good (x^ < 2.0) and the residuals 
are noise-dominated. This supports the conclusion that galaxies 
without significant structure result in acceptable fits irrespective 
of whether they are binned or not. 



vestigate correlations between galaxy colours, effective sizes 
and concentrations. The (inverse) concentration index was 
found to correlate well with the visual estimates of morphol- 
ogy. Shimasaku et al. (2001) have kindly provided us with 
their visual morphologies in order to compare them with the 
Galactica S/T ratios. The are 166 galaxies in common in 
the two samples which have xt < 2.0. Fig. IC3I shows that 
there is a fair correlation between the (inverse) concentration 
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Figure B6. The distribution of R50 (the radius enclosing 50% 
of the Petrosian flux) for galaxies in the SDSS which match our 
selection criteria (solid histogram). Results are also shown for the 
sample used by Tasca Sz White (2006) (dotted histogram). For 
comparison, we indicate the maximum allowable seeing for our 
sample by the vertical dashed line. 




Figure CI. Galactica S/T ratio vs. {u—r) colour for the sample 
of 7493 SDSS galaxies studied here. Most galaxies colour-classifled 
as late types are predominately disk-dominated systems while 
most galaxies colour-classifled as early type have S/T > 0.4. 



index and the visual morphology for these 166 galaxies, con- 
firming the conclusions of Shimasaku et al. (2001). Plotted 
in Fig. IC4l is the correlation between our derived S/T ratios 
and the visual morphology for these galaxies. Although the 
scatter is large, there is a clear trend for the earlier types to 
have larger S/T ratios. 
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Figure C2. Galactica S/T ratio vs. (inverse) concentration in- 
dex, C = R^o/Rqo, for our sample of 7493 SDSS galaxies. More 
centrally concentrated galaxies (C < 0.33) predominately have 
higher S/T ratios. 




Figure C4. Galactica S/T ratio vs. the visual morphology de- 
termined by Shimasaku et al. (2001) for the 166 galaxies in com- 
mon in the two sample. There is a general trend for the S/T ratio 
to increase along the S-SO-E morphological sequence but the scat- 
ter is large. 
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Figure C3. Concentration index vs. eye morphology for 166 
galaxies in common between our sample and that of Shimasaku 
et al. (2001). 



